Role of Divergence of Classical Trajectories in Quantum Chaos 
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We study logarithmical in h effects in the statistical description of quantum chaos. We found 
analytical expressions for the deviations from the universality in the weak localization correction and 
in the level statistics and showed that the characteristic scale for these deviations is the Ehrenfest 
time tE = A -1 1 ln7i|, where A is the Lyapunov exponent of the classical motion. 
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It is accepted in literature to call the consideration 
of quantum phenomena in classically chaotic system the 
"quantum chaos" Q. For the de Broglie wavelength A_f 
much smaller than the characteristic size of the system, 
the quantum phenomena still bear essential features of 
the classical chaotic motion. Examples of such systems 
studied both theoretically and experimentally are ballis- 
tic cavities or antidot arrays The quantities usually 
considered include different correlators of quantum spec- 
tra of the system (level statistics) as well as of different 
response functions, e.g. fluctuations of the conductance 
(mesoscopics) or the quantum correction to the averaged 
transport coefficients (weak localization). 

In principle, all the aforementioned characteristics can 
be found by solving the one particle Schrodinger equation 
for the given system. However, the Schrodinger equation 
for such systems can not be solved analytically. Substan- 
tial progress can be achieved in the statistical approach 
to the quantum chaos. In such approach one gives up 
attempts to find a contribution of a single quantum state 
but instead studies correlators averaged over large num- 
ber of quantum states. The averaging for a given system 
can be performed either over wide range of energies or 
over the applied magnetic field. 

In the present Communication we apply the supersym- 
metry description H || to investigate how the univer- 
sality is established in the statistical properties of the 
system at low frequencies. We will show that the time 
it takes to establish the universality is 4e = A _1 |ln7i|, 
where A is the Lyapunov exponent of the classical mo- 
tion. We will express deviations from the universality 
in the level statistics and in the weak localization cor- 
rections in terms of the single renormalization function 
||. Finally, we emphasize the necessity of the finite re- 
qularizator in the Perron- Frobenius operator in obtaining 
physical results for the quantum corrections. 

Let us first discuss the physical origin of the logarith- 
mic in h corrections. In the semiclassical approximation, 
each classical trajectory corresponds to the quantum me- 
chanical amplitude. Quantum phenomena in the system 
originate from the interference of the different ampli- 
tudes. After the averaging, the most of the interference 
contributions vanish. Only ones which survive are the 



products which contains the pairs of the coherent ampli- 
tudes. Such coherent amplitudes are contributed by the 
segments of the same classical trajectories and the result- 
ing products are expressed in terms of the probabilities 
to find such segments. The latter probabilities are found 
by solving the classical equation of motion. Most usable 
quantities are the probability where the initial i and final 
/ states coincide nj = r^, Yj = [we will denote this 
probability as X>+(t;nj, r^)] or related to each other by 
time inversion rif — — n^, r/ = [we will denote this 
probability as T>_(t; iij, r$)]. Here r and n are the coor- 
dinate of the particle and the direction of its momentum 
respectively. First quantity is relevant for the leading 
approximation for the two point correlator of the density 
of states (DoS) [||J^J|,§1, whereas, the second quantity 
is important for the weak localization correction to the 
conductivity Q and for the higher order approximations 
for the correlator of DoS, see below. 

In what follows we will consider only ergodic systems. 
This means that after some time the particle visits all 
the phase space allowed by the energy conservation, i.e. 
the classical probabilities T>± averaged over the condi- 
tions rii, Yi cease to depend on time and take value of 
1/S, where S is the volume of the system. It is very cru- 
cial that the equilibration time for 2?_ is parametrically 
larger than that for the probability T>+ . 

Characteristic relaxation time for the probability 2?+ 
is of the order of the flying time of a particle across the 
system r/j ~ L/vf for the ballistic regime or the Thou- 
less time tt — L 2 /D for the diffusive regime, (L is size 
of the system, D is the diffusion coefficient, and Vf is the 
Fermi velocity). 

On the other hand, in the strictly classical limit prob- 
ability to have n/ = — n;, Tf = vanishes no matter 
how large traveling time t is. This is due to the fact that 
the final state can be reached by moving along a classical 
trajectory which coincides with the initial one. It means 
that a particle must be reflected exactly backwards from 
an obstacle. For the chaotic system, the measure for 
such process is zero and that is why 2?_(t;rij, r^) = 0. 
The only reason for this probability not to vanish is 
that the initial and final conditions can not be speci- 
fied with the accuracy better than it is allowed by the 
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uncertainty principle. Due to this principle, the dif- 
ference |n/ x n,j| = 6(f>o can not be smaller than the 
diffraction spreading 8<fro > y/Xp/a, with a 3> Xp be- 
ing the characteristic spatial scale of the static potential 
the particle moves in. In order to find the probability 
for such close (but not coinciding) r/, r^, one has to 
take into account the fact that the motion of the par- 
ticle at the initial and final stages are correlated. This 
is because the trajectory along which the particle moves 
on the final stage, [r(i — ti), ~n(t — ii)] almost coincides 
with the trajectory particle moved along at the initial 
stage, [r(ti),n(£i)]. This problem is equivalent to the 
consideration of the divergence of two classical trajec- 
tories ("1" and "2") which start from the same point r$ 
with a small difference in the directions of their momenta 
1 112(0) x ni(0)| = Sfio, (it can be seen by the time inver- 
sion on the final segment). In the chaotic system the dif- 
ference <50i2 (£) = |ri2(i) x ni(f)| grows exponentially with 
time 8(j)(t) ~ 8<f>Qe xt where A is the Lyapunov exponent of 
the classical chaotic motion. Therefore, we have also for 
the given trajectory 8<j)(t\) = |n(t— 1\) Xn(ti)| ~ 8cj) e xtl . 
In order to close the trajectory at some time t\ < t/2, 
angle 8(j)(t\) should become of the order of unity and thus 
t > (2/A) In (l/Scf)o). Taking into account 8(j)o ^ \JXp/a 
we conclude that the time it takes to establish the equi- 
librium value of function 2?_ is the Ehrcnfcst time 
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and T>- vanishes at smaller time, X>_ ~ (t — tp). 

The above discussion leads us to the following expres- 
sions for the Fourier transform of the classical probabil- 
ities T>± for the frequencies uj smaller than inverse time 
of the travel of the particle across the system 
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where uj + = u> + iO. The denominators in Eqs. (|2|) reflects 
the ergodicity of the system at large time and the renor- 
malization function T(u>) describes the delay of T>_(t) 
with respect to T>+(t) by the Ehrenfest time tp 



T(u) — exp iujtp 



ui 2 X 2 t E 
X 2 



(3) 



The second factor in Eq. (j^) characterizes the fluctua- 
tions of the Lyapunov exponent and the parameter A2 is 
of the order of A. More details about the derivation of 
function T can be found in Ref . M . Appearance of the 
new time scale tp is the qualitative difference between 
quantum chaos a > Aj? and quantum disorder a < Xf 
regimes. (In the systems where the scale of the potential 
a is different from the transport mean free path l tr , the 
criteria for quantum chaos is a ^> Vh^Xp, see Ref. M ) . 



A powerful method for the calculation of the averaged 
quantities is the supersymmetric nonlinear er-model pi- 
oneered by Efetov |3) for the disordered systems. Re- 
cently, the supersymmetric action was suggested by by 
Muzykantskii and Khmelnitskii Q and more recently by 
Andreev et. al. I| for the system in the ballistic regime. 
Effective action in Ref. || is defined by means of the clas- 
sical Perron-Frobenius operator which differs from the 
first order Liouville operator by the regularizator of the 
second order. This approach enables one to perform the 
systematic semiclassical expansion for the averaged quan- 
tities and to understand how the underlying classical dy- 
namics shows up in various quantum correction. 

Partition function Z in the supersymmetry approach 
is given by the functional integral || 



Z{J}= / DQ(l)exp 



~ J dlSTr(£ + £j) 



(4) 



with Lagrangians C, Cj being defined as 



= T AT 



Cj = »Q(1) 



./, A A f ( J 2 (I)A+ + J 2 *±l(l)A. 
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where v is the density of states per unit area. We used the 
short-hand notation 1 = (ni,ri), 1 = (— ni,ri), dl = 
dniG?ri/27r, coordinate r and the direction of the mo- 
mentum n = (cos 4>, sin (f>) characterize the position of 
the particle on the energy shell, (we will restrict our- 
selves to two dimensional systems) . Liouvillean operator 
L describes the classical evolution on the energy shell and 
defined by the Poisson bracket £• = {•, H}, where Ti. is 
the Hamiltonian function. Last term in the Lagrangian 
C is the regularizator physical significance of which will 
be discussed later. Operation of supertrace is defined in 
Ref. H . We will consider only systems with the unbroken 
timereversal symmetry (orthogonal ensemble). 

In Eq. (|^), f is 8 x 8 supermatrix defined in a linear 
superspace p ® g ® d which we represent as the direct 
product of three linear spaces; p and d are the spaces 
of retarded-advanced and time reversal (complex con- 
jugate) 2x2 matrices respectively, and g is the super- 
space of fermion-boson 2x2 supermatrices. All the rel- 
evant matrices can be conveniently expressed in terms 
of the Pauli matrices r", = (r° ± zt^)/2, acting 
in spaces a — p,g,d. Matrices A = rf ® l 9 ® l d , 



I 9 ® l d and k = V 



l d break the 



symmetry in advanced-retarded and in fermion-boson 
spaces respectively. Matrix T is the subject to constraints 
T^KT = K and T\l) = CT T (l)C T 7 where 2K = 
[(1 + rf) <g) I 9 + (I - rf ) ® r|] <g) l d and matrix of charge 
conjugate C is given by C = l p <g> (l ff (g) — t% (g) 7+). 

Partition function (|J) allows to find different quan- 
tum mechanical correlators averaged over the wide range 
of energies. Two point correlator of the DoS, R(w) — 
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A 2 (p(e + uj)p(e)) e (where p(e) = Tr5(e-H), with H be- 
ing the Hamiltonian of the system, and A is the mean 
level spacing ) and two-particle Green function T> can be 
found as certain derivatives of action 



2?( W ;1,2) = -^ 



A 2 9*Z 



16tt 2 9 2 Ji 
S 2 Z 



v SJ 2 (l)SJ 2 (2) 



Jl,2=0 



Jl,2=0 



(6a) 
(6b) 



We will be interested in effects associated with the en- 
ergy scale u) ~ t^ 1 3> A, and will not consider contri- 
butions oscillating on A. For analysis in this situation it 
suffices to consider only small fluctuations of Q around 
A. We use the standard parametrization 



T= 



1+iP 

7T+F' 



PA = -AP, KP^K=P, P(l)=-P(l), (7) 



where the operation of time reversal is defined as 

M = KCM T C T K, (8) 
for arbitrary supermatrix M. 



Substituting formula Eq. (fjj) into Eq. (Eh and keeping 
terms up to the fourth order in P we obtain C = £o+A'n* 
where the quadratic part of the Lagrangian describing the 
classical motion is given by 

C a = P(-iuj + L R )P, L r = L--^-. (9) 



Operator L R is known as the Perron-Frobenius opera- 
tor. Quartic part, responsible for lowest order quantum 
corrections has the form 



1 / dP 

C mt = -P 3 {-iu + L R )P + - (P— 



(10) 



We expand the partition function up to the first or- 
der in Ci n t and up to the second order in the source 
Lagrangian 

Cj = iJ x k (1 - 2P 2 + 2P 4 ) + 2J 2 (k + 1) A+ (P - P 3 ) , 

where we used Eq. (R) and omitted traceless terms. For 
calculating the averages of the arising products of ma- 
trices P we use the Wick theorem with the contraction 
rules 



2ttvP (\)M P{2) =V°(l,2)AfSTr 



MA 



■V°(l,2)Au STr 



MA+ 



X>°(l,2)A | j"MA+ +P°(l,2)AjJ"MA | 7, (11 



27nASTr[MP(l)] STr[NP{2)} = STr (V°(l, 2)M — V°(l, 2)M) A^NAj + (T>°(1, 2)M — 2)M) At A^Ajj 



where matrices A|J : = (l±A)/2 break the symmetry in 
the retarded-advanced subspace, and the classical prop- 
agator T>° is the solution of the equation 



(-iuj+ + L B }j ^ 



-iuj+ + L r ) X>° (1,2) = 27r*ia. 



(12) 



Substituting the result of the averaging in Eq. (pb|), we 
obtain P(l, 2) = 2) + 5T>(1, 2), where the quantum 

correction to the classical propagator is given by 

2nv5V(l, 2) = D°(l, 2)V°(2, 2) + D°(l, 1)P°(1, 2) (13) 
dSD°(l, 3)X>°(3, 2) (2icv - L R ^V a {i, 3). 



Notice, that Eq. (13) gives a correction only to 
the non-zero modes of the Perron-Frobenius operator, 
J dl#D(l,2) — 0, which is a consequence of the charge 
conservation. 



Analogously, we obtain with the help of Eq. (6a) the 
following expression for the two-point correlator of DoS, 
R = R° + SR. Here 



A 2 r 

R°(uj) = 1 + — Im / dl 



dV a {l,l) 



(14a) 



which is well-known result for the disordered and 
chaotic systems Quantum correction, obtained for 
the first time in the present paper, has the form 



5R(u) 



2n 3 . 



-Re 



4tw Jdld2V°(l 1 2)5V°(2,l) 



r/3 



92?°(3,3) 
dui 



2iuj — L 



<9P°(3,3) 



(14b) 



acts on both ar- 
we used identity 



Operator L R in Eqs. ( |13| ) and (J14F 
guments 3,3. Deriving Eqs. (|14|). 
-id u V°(l,2) =J d3D°(l,3)D (3,2).' 

Equations ( p^| ) and ( |l4b| ) describe the lowest quan- 
tum corrections expressed in terms of the solutions of 
the Liouville equation (with the regularizator added) for 
a given system, where no ensemble averaging is assumed. 
They are both determined by the classical propagators 
between the points related by the time inversion. As we 
already explained, the equilibration time for such proba- 
bilities is the parametrically large Ehrenfest time. 
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Let us first discuss the two-point DoS correlator R. 



The classical propagator T>° entering into Eq. (14a) is 
the probability of the return to the initial state and at 
energies uj much smaller than the Thouless energy it co- 
incides with V + from Eq. (||). First term in Eq. ( |l4b| ) 
correspond to the weak localization corrections of non- 
zero modes of Perron-Frobenius operator and it can be 
neglected at such low-frequencies. Propagators entering 
into the second term are the classical probabilities with 
the initial and final states related by the time inversion 
and they coincide with 2?_ at low frequencies. Substitut- 
ing Eqs. H) into Eqs. (113), we obtain 



R(uj) 



1 



A 3 w ( d r(w) 



-Im 



(15) 



where the renormalization function T defined in Eq. (|^). 
It is easy to see that result (|l5]) does not contain terms 
linear in tE- Actually, this can be proven for all orders 
of perturbation theory in A/lo by using the approach 
similar to Ref. H). Even though the found correction 
to the universal Dyson result |hJ for the orthogonal en- 
semble is small, it oscillates with the period t E , where 
A <C t E x <C Tji and, therefore, can be distinguished. 

The result for the quantum correction to the conduc- 
tivity of infinite system 8a obtained from Eq. ( |l3| ) 

a2 ( 



In 



1 



2-K 2 h \ljJTtr 



r*( w ). 



is renormalized in comparison with the quantum disor- 
der regime |l^] by T. Details of the derivation of this 
equation from Eq. (13) can be found in Ref. §. 

Ehrenfest time (|1) contains de Broglie wavelength Xp. 
This scale is already absent in the effective er-model (||) 
which is formulated on the Hilbert space of functions 
smooth on scale A f , and the lower cut-off of the logarithm 
is related to the reqularizator 1/r in the Lagrangian (|^). 
It follows from the solution of Eq. (|l2|), [g that P_ has 
the form (J2j) with tp = A -1 In (At). Therefore, the reg- 
ularizator can not be put to zero even in the end of the 
calculation and it should be assigned some physical value. 
The value of the reqularizator can be derived for the case 
if, in addition to the semiclassical potential, there are 
also quantum impurities in the system which provide the 
small angle scattering. The diffraction on the semiclassi- 
cal potential itself is described by an equation more com- 
plicated than (|lj) and the regularizator was not found 
consistently within the a-model approach. However, it is 
not really necessary because the dependence on the regu- 
larizator is only logarithmical. Using Eq. (|j) obtained by 
different arguments, we conclude that the physical value 
of the regularizator is given by 



A—. 

a 



(16) 



There is a subtlety in Eq. (|13|) which deserves more dis- 
cussion because it helps to understand the importance of 



the reqularizator in the supersymmetric Lagrangian (|^). 
With the help of Eq. (|l|), one can rewrite Eq. ( |l3] ) to 
the more compact form 



SV(1,2) = J 



(IS 



P°(3,3) <9P°(1,3) <9P°(2,3) 



Naively, 8T> — > for r — > oo. However, this contribution 
is anomalous and caution should be exercised while tak- 
ing the limits. Namely, D°(l, 3 + S(/))V (2, fa - 5(f)) is a 
singular function on 5(j) and this singularity is cut-off by 
the same regularizator 1/r. As the result, the derivatives 
over cj) s are proportional to yfr and the dependence on the 
regularizator remains only logarithmical, see Eq. (fi3|). 

In conclusion, we studied logarithmical in h effects in 
the statistical description of quantum chaos. We found 
analytical expressions for the deviations from the univer- 
sality and showed that the characteristic scale for these 
deviations is the Ehrenfest time tp = A _1 | lnfi.|, where A 
is the Lyapunov exponent of the classical motion. Finally, 
we discussed role of anomalies in the supersymmetric a- 
model ||] °f quantum chaos. 
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